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Preconditioning Maximal Center Gauge with Stout Link Smearing in SU(3) 
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Center vortices are studied in SU(3) gauge theory using Maximal Center Gauge (MCG) fixing. 
Stout link smearing and over-improved stout link smearing are used to construct a preconditioning 
gauge field transformation, applied to the original gauge field before fixing to MCG. We find that 
preconditioning successfully achieves higher gauge fixing maxima. We observe a reduction in the 
number of identified vortices when preconditioning is used, and also a reduction in the vortex-only 
string tension. 

PACS numbers: 12.38.Gc 11.15.Ha 12.38.Aw 



I. INTRODUCTION 

Despite more than 30 years of intense study, quark 
color confinement in hadron physics remains unexplained 
(for a recent overview see Ref. [l|). Within the frame- 
work of lattice gauge theory, the prevailing view is that 
quark confinement is the result of a particular class of 
gauge field configurations which dominate the QCD vac- 
uum on large distance scales. Two potential candidates 
have been most commonly investigated: confinement by 
means of Zn center vortices and confinement due to 
Abelian monopoles (for a critical discussion of both sec 
Ref. 0). To enhance these particular features, gauge 
fields can be first fixed to a suitable gauge, such as 
Maximal Abelian Gauge (MAG) 3] or Maximal Center 
Gauge (MCG) [1]. Monopoles and center vortices are 
then defined by the projection of these gauge-fixed fields 
onto /7(l) Ar_1 or Zn, respectively. Significant progress 
to date has occurred in SU{2) using MAG and MCG, 
with original findings reproducing about 90% [5] and 
about 100% [H, respectively, of the non- Abelian string 
tension. Remo ving monopole @, [H, [l3| or center-vortex 
d 03, El El El El degrees of freedom from SU(2) 
lattice gauge fields appears to leave topologically triv- 
ial, non-confining gauge fields that do not spontaneously 
break chiral symmetry. 

The significance of the center of the gauge group is 
what connects possible candidates for this special class 
of configurations. As outlined for the case of Lapla- 
cian Center Gauge (LCG) in Ref. [15], all monopole 
world lines are embedded in 2-dimensional vortex sur- 
faces. These topological objects naturally occur together 
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as local gauge defects. In MCG it has been observed 
computationally that over 90% of monopole currents are 
localized on center vortices 

[ll, E3- Strongly correlated 
effects between the two have also been observed by means 
of studying monopoles after vortex removal and vice- 
versa [18|, as well as through the effect of their removal 
on the spectra of the overlap Dirac operator (l3| . 

Again, all these advancements have been in SU(2) 
and work in SU(3) has not progressed to this level. 
While initial investigations were hopeful [HI, Ho|, sub- 
sequent results for MAG and MCG were not so encour- 
aging [U [H, [23|, as they failed to reproduce the full 
non-Abelian string tension. Earlier studies in SU(2) us- 
ing MCG reported that the center-projected configura- 
tions recovered the full string-tension, however further 
study into the ambiguities of the gauge-fixing procedure 
showed that this result is plagued by Gribov copy effects 
[24l I25I [26j : methods which give higher values of the 
gauge fixing functional produce smaller values for the 
vortex-induced string tension. We point out that when 
the Laplacian Center Gauge of Refs. [HI, [27j (which is 
free of Gribov ambiguities on the lattice) is used as the 
fixing method, the full SU(3) (and SU(2)) string tension 
is recovered for the center-projected gauge fields but only 
in the continuum limit. However, unlike MCG vortices 
(2c| , the interpretation of LCG vortex matter is cumber- 
some in the same continuum limit [22L I29I ] . 

In this paper we focus on the Gribov problem of the 
SU(3) center-vortex picture of confinement using the 
MCG fixing method. We apply the "smeared gauge 
fixing" method of Ref. [30] to MCG to ameliorate this 
Gribov problem. This creates a pre-conditioning gauge 
transformation for the configuration that should bring it 
closer to the global maximum. We investigate the effect 
of this method on the features of the long-distance behav- 
ior of the static quark potential as evaluated on configura- 
tions where the P- vortices derived from MCG have been 
removed and configurations composed purely of these 
P-vortices. In SU(2), it has been shown that center- 
vortex removal specifically targets topological properties 
[Til EH , so as well as using stout-link smearing j3l| we 
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also employ over-improved stout-link smearing [32j to at- 
tempt to exploit the link to topological structure (33j . 



II. METHODOLOGY 
A. Identifying Vortex Matter 

In the center-vortex picture of confinement the gauge 
fields are considered to be decomposed into a long-range, 
smooth field Z^ carrying all the confining fluctuations 
and a short-range field Vu containing non-confining per- 
turbations as well as other short-range effects 

Ufa) = Zfa)Vfa). 

Here Z^ix) is the center element which is closest, on the 
SU(3) group manifold, to U^{x). A vortex is a configu- 
ration of the gauge potentials topologically characterized 
by non-trivial elements of Z% and is created by a singu- 
lar gauge transformation. The non-trivial center element 
of the singular gauge-transformation characterizing the 
vortex may be made to be distributed over many links 
of an encircling loop (due to the short-range effects of 
V fl (x)). If we assume that by a gauge transformation the 
non-trivial center element can be concentrated on just 
one link we can compress this thick vortex into a thin 
one. If we then project this gauge transformed config- 
uration onto its center elements, the projected vortices 
(P-vortices) linking with the loop should then correspond 
to the thin vortex. It is for this reason that we adopt the 
use of gauge-fixing to obtain the necessary gauge trans- 
formation. It is the choice of gauge that determines our 
method for finding the center vortices and, therefore, the 
connection between the P-vortices and the thick center 
vortices present in the original configuration. The par- 
ticular choice of gauge, the properties of the P-vortices 
associated with each choice and the Gribov problem that 
it creates is what has polarized opinions in this area 

EE in m, m, H3i . 

Here, we employ the MCG gauge-fixing algorithm as 
outlined in Ref. [22]. The gauge condition we chose to 
maximize (with respect to the center elements Z ll (x)) in 
this algorithm is 

where Ni is the number of links on the lattice and U n is 
the gauge-transformed field. 

After fixing the gauge, each link should be close to a 



center element of SU(3), Z 



1 



€ { — 1,0, 1}. Since, for every link, 



trC/"(x) = u^e^"-" and <f> Xtll = tan 



= -^-m with 



_ 1 ImfoCff (*)) 
Re(trL^(ai)) 



then <fi x „ should be near to some <j) m , by construction of 
the gauge-fixing condition. We then perform the center 



projection by mapping 

SU(N) i-> Z N : Uf L {x) h-» Zfa) with Zfa) = e**-", 

with the appropriate choice of <f>™„, m G { — 1, 0, 1}. 

To reveal the vortex matter we simply take a product 
of links around an elementary plaquette. We say a vortex 
pierces the plaquette if this product is a non-trivial center 
element and the plaquette is then a P-vortex. We can 
remove these P-vortices by hand from the configuration 
using Ufa) = Zfa)Uj}(x). 



B. Smearing as a Preconditioner 

In the center- vortex picture of confinement, the center 
elements correspond to the long-range physics. It would 
seem reasonable then to employ the use of smearing to 
smooth out the short-range fluctuations and allow the 
gauge transformation to see more of the underlying long- 
range physics. To this end we construct a precondition- 
ing gauge transformation for each gauge field to obtain 
higher maxima in the gauge-fixing procedure and thereby 
directly address the Gribov-copy issue [30l ]. 

Firstly, we smear the gauge field using any smearing 
algorithm (stout-link smearing (3l| has been applied here 
as well as over-improved stout-link smearing which has 
been shown to better preserve the topological structure 
underlying the original configuration [32]). We then fix 
the smeared field using the MCG gauge-fixing method. 
At each iteration we keep track of the total gauge trans- 
formation that has been applied to the smeared gauge 
field. Once the algorithm has converged we use the stored 
total transformation as a preconditioning gauge transfor- 
mation for the unsmeared gauge field. We emphasize that 
the (unsmeared) preconditioned gauge field remains on 
the same gauge orbit since the preconditioning is merely 
a (specific) gauge transformation on the original links. 
Gauge-fixing the preconditioned field simply gives us a 
Gribov-copy of the result from gauge-fixing the original 
gauge field. 



III. RESULTS 

Calculations are performed using 200 quenched config- 
urations with the Liischer-Weisz plaquette plus rectangle 
gauge action [3f| on a 20 3 x40 lattice with (3 = 4.52. Sim- 
ilar preliminary results have being found on 100 16 3 x 32 
lattices (with (5 = 4.6) and have been reported elsewhere 

Stout-link smearing with a smearing parameter of 0.1 
is used to construct the preconditioning transformation 
with the number of sweeps ranging from to 20 in steps 
of 4 sweeps. We also employ over-improved stout-link 
smearing with a smearing parameter of 0.06 and an e pa- 
rameter of —0.25. Here, each preconditioning was con- 
ducted independently. 
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TABLE I: 


Results for the average maximum gauge condition 


16 ±9 


14±8 



Vt/[0] as a function of preconditioning stout- link smearing 
sweeps (OI signifies over- improved stout-link smearing). For 
each of the sweeps used in the preconditioning: the average 
total (smeared gauge field fixing plus preconditioned gauge 
field fixing) number of blocks of 50 iterations used, the aver- 
age smeared gauge condition maximum reached, the average 
preconditioned gauge condition maximum reached and the 
percentage of plaquettes that are P-vortices. 



Given that the original goal was to increase the gauge- 
fixing maxima achieved in MCG fixing, we can see from 
Table U that we are successful, in this regard, in every 
case. With each level of preconditioning a higher gauge 
condition maximum is achieved both for the smeared 
gauge field and the preconditioned original field. If we 
compare and 4 sweeps of preconditioning, we can see 
that the magnitude of this increase is initially large but 
the increase is slower as we precondition to higher lev- 
els. However this increase does not come without a cost, 
the number of gauge-fixing iteration blocks (a block is 
50 iterations) required almost doubles between the un- 
preconditioned fixing and the maximum amount of pre- 
conditioning. Typically, two-thirds of the iterations are 
spent fixing the smeared field and one-third fixing the 
preconditioned field. 

What is most significant about this table however is 
that with each level of preconditioning the percentage of 
projected plaquettes which are P-vortices drops signifi- 
cantly. Without preconditioning 3.21% of plaquettes are 
vortices and this drops to as low as 1.28% for the highest 
level of preconditioning. 

In Table [TT] we investigate further by looking at this 
particular effect between all the different levels of precon- 
ditioning. In the upper triangle of this table (from pre- 
conditioning level row to preconditioning level column) 
we report the percentage of configurations that experi- 
ence a reduction in the measured number of P-vortices. 
As we can see, this percentage is always high but the 
effect is lessened as we move to transitions, particularly 
small transitions, between higher levels of precondition- 
ing. It should be noted however that the relative differ- 
ence between, for example, 20 sweeps of stout-link and 
80 sweeps of over-improved stout-link preconditioning is 
difficult to quantify but the effect is still significant for 
this transition. 

The magnitude of this effect is also reported in Ta- 



TABLE II: Comparisons between different preconditioning 
levels of stout-link smearing (OI signifies over-improved stout- 
link smearing). The upper triangle of this table (from precon- 
ditioning level row to preconditioning level column) we report 
the percentage of configurations that experience a reduction 
in the measured number of P-vortices. The lower triangle 
(preconditioning level column to preconditioning level row) 
of this table gives the percentage reduction of the number of 
P-vortices for the configurations that experienced a reduction. 



ble[TTJ When reading from the lower triangle (precondi- 
tioning level column to preconditioning level row) of this 
table we can see the percentage reduction of the num- 
ber of P-vortices for the configurations that experienced 
a reduction. In the transition from no preconditioning 
to any other level, the order of a 50% reduction is ob- 
served. For other transitions it would seem the effect 
drops to the 10% level reasonably quickly, but again we 
see an increased effect when we consider over-improved 
smearing. It should be noted that, regardless of the pre- 
conditioning level, the center phases of the links of the 
fields always remain evenly distributed across the three 
possible values, reflecting the fact that the realization of 
center symmetry remains unaffected. 

We can look to Table IIIII when considering whether 
a higher gauge-fixing maximum translates into a lower 
number of P-vortices. When reading from precondi- 
tioning level row to preconditioning level column, the 
percentage of configurations that experience an increase 
in the gauge-fixing maximum is shown. Similar trends 
to that of Table [Til ar e observed, with large effects ini- 
tially which become reduced for small transitions be- 
tween higher levels. Of these configurations with an in- 
creased maximum we can see almost exclusively (when 
reading from preconditioning level column to precondi- 
tioning level row) that an increased gauge-fixing maxi- 
mum does lead to a lower number of P-vortices. 



A. The Static Quark Anti-quark Potential 

The fact that we can reduce the number of P-vortices 
through preconditioning is not necessarily a cause for 
concern. As stated previously, our method for determin- 
ing the location of center vortices is only justified by the 
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FIG. 1: The effective potential plots for the lowest (left) and highest (right) levels of preconditioning smearing. The upper plots 
contain the data for the original gauge-fixed configurations, the middle plots contain the data for the vortex-only configurations 
and the lower plots contain the data for the vortex-removed configurations. Each plot contains data for a range of quark 
separations. 
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TABLE III: Comparisons between different preconditioning 
levels of stout-link smearing (OI signifies over-improved stout- 
link smearing). When reading from sweep row to sweep col- 
umn (upper triangle) the value shown is the percentage of 
configurations that achieve a higher gauge fixing maximum. 
When reading from sweep column to sweep row (lower tri- 
angle) the value shown is the percentage of the configura- 
tions with higher maximum that achieve a lower number of 
P-vortices. 

physical relevance of the P-vortices that we determine. 
A first step in determining this relevance is the calcu- 
lation of the static quark anti-quark potential. In the 
center- vortex picture, the string tension a as determined 
from the infrared behavior of this potential should be 
fully accounted for by the center-vortex component of 
the gauge fields, Z^, with the Coulombic term accounted 
for by the vortex removed component, V^. Since we 
can "remove" the determined P-vortices by the operation 
U' (x) = Z^(x)Uj}(x), we can seek to observe these prop- 
erties directly. However, since the determined P-vortices 
are gauge-dependent (and their number Gribov-copy de- 
pendent, as we have already seen) then so too are the 
subsequent measurements of the static quark potential 
from the vortex-only and vortex-removed components of 
the configuration. 

Computing the static quark anti-quark potential as a 
function of the quark separation is a two step process. 
Wilson loops W(R : T) of extension of R x T have the 
large T behavior 

(W(R,T)) (xexp{-V(r)aT}, r := Ra, 

where a is the lattice spacing. The method for extracting 
the effective potential is thus identical to that of extract- 
ing effective masses from two-point functions in hadronic 
spectroscopy. To obtain the static quark anti-quark po- 
tential as a function of the quark separation we simply 
repeat this process for a range of values of the separa- 
tion R. By using off-axis spatial paths for the Wilson 
loops, we can obtain non-integer values of R. We exploit 
full space-time translation to improve the statistics of our 
Wilson loops. 

Since the final plot is composed of fits performed on a 
large number of effective potential plots for all the differ- 
ent separations, it is prudent and necessary that the fac- 



tors determining those fits are given, and taken into ac- 
count, when analyzing the subsequent static quark anti- 
quark potential as a function of separation. The difficul- 
ties associated with such fits can be easily recognized in 
Fig. [TJ In these plots we show the static quark potential 
for a variety of quark separations for each of the original, 
vortex-only and vortex-removed configurations. On the 
left we show these plots for the unpreconditioned MCG 
fixing and on the right we show the same plots for 80 
sweeps of over-improved smearing as a preconditioner. 

One of the first things to discuss is the difficulty in ob- 
taining a satisfactory fit range for the data, particularly 
in the case of the vortex-removed configurations. With 
these configurations, more so at larger separations, the 
potential falls rapidly and decays into noise quickly. A 
visually satisfying plateau region is not evident and we 
must rely on the fitting routine to determine the good- 
ness of the fit. What the plot can tell us is that the 
effective potential continues to fall (at separations of 5 
lattice spacings and greater) until at least time-slice 5. 
Since the data decays into noise around this point, we 
chose to constrain our fit using timcslicc 5 and fit from 
this slice to slice 7 (a straight-line fit to 3 points). This 
constraint is then applied to all values of the separation. 
What we find is that while this may lead to reasonable 
(< 1.3) values of the \ 2 P er degree of freedom in the ma- 
jority of cases, there are certainly significant deviations 
from this desirable result. 

Global fit ranges are chosen in a somewhat similar way 
for the unpreconditioned and vortex-only configurations. 
For the unpreconditioned configurations, the global fit- 
range was chosen to be between timeslice 4 and 6 since 
these accounted for the systematic drift of the potential 
at large separations for lower time values and also gave 
reasonable \ 2 behavior. In the case of the vortex-only 
configurations, the errors are far more controlled but the 
potential rises at small times and plateaus far later so 
the fit range was chosen to be from timeslice 10 to 12 
but again some of the \ 2 P er degree of freedom values 
were unsatisfactory. This is most likely due to the heavy 
constraints placed on the fit by the accurate potential 
determinations. 

Of significant concern when comparing the potentials 
of the unpreconditioned and preconditioned results in 
Fig. [T] is the direct comparison of the potentials in each 
case. The top figures show the potential after gauge- 
fixing but prior to the center projection and vortex- 
removal and, of course, since the static quark potential 
is gauge invariant these plots are identical. For the low- 
est plots, which contain the vortex-removed data, we can 
see that the rate of decay for the preconditioned data 
has dropped and the quality of the data does not al- 
low us to see whether it is possible that it plateaus at 
the same level as the unpreconditioned data. Of course 
the most dramatic effect occurs in the middle plots with 
the vortex-only data. There is a dramatic reduction in 
the magnitude of potential for all separations and this 
is direct manifestation of the Gribov-copy effect for this 
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FIG. 2: Static quark anti-quark potential plots for each of preconditioning smearing sweeps used. Each plot contains data for 
the full, vortex-removed and vortex-only configurations. The data shown uses a 3 timeslice fit window in each case with the fit 
window being from timeslice 4 to 6 for the full data, timeslice 5 to 7 for the vortex-removed data and timeslice 10 to 12 for the 
vortex-only data. 
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FIG. 3: The static quark anti-quark potential plots for both the lowest (left) and highest (right) levels of preconditioning. Only 
the on-axis data is shown for the vortex-removed configurations. 
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FIG. 4: Ratio between the vortex-only string tension and the 
vortex density as a function of preconditioning. 

gauge-fixing method. 

This Gribov-copy effect is also manifest in Fig. [51 Here 
we show plots of the effective potential as a function of 
separation for each of the six levels of preconditioning 
used as well as the unpreconditioncd data. It would ap- 
pear that the findings are consistent with loss of confine- 
ment upon P-vortex removal. Although it would seem 
that this is perhaps not such a reasonable observation in 
the over-improved case, if we look exclusively at the on- 
axis contributions to the potential in this case (Fig. [3]) 
and compare it to that of the unpreconditioned case, we 
observe that a plateau in the potential may exist but at 
larger values of the separation. This would concur with 
our previous observation that the potential takes longer 
to plateau in this case and therefore the fit window may 
not be adequately account for this effect. Careful ex- 
amination of the vortex-removed plot also reveals that 
we obtain an increasingly more accurate fit to the short- 
range Coulombic portion of the potential. 

What is more significant however is that the value of 
the string-tension determined from the vortex-only con- 
figurations drops dramatically, and systematically, from 
~ 60% to as low as ~ 16% of the full string tension with 
increased preconditioning providing improvement. This 
is a disturbing manifestation of the Gribov problem since 
it perhaps questions how accurately we have determined 
the center vortices by our projection of the P-vortices 
with our fixing method. 

B. Discussion 

The use of smearing as a preconditioning technique 
does indeed lead to higher maxima in the MCG gauge- 
fixing condition Vjy[f2]. These higher maxima in turn 
lead to lower numbers of P-vortices determined in the 



center projection. In SU(2), similar results have been 
obtained when seeking higher maxima by means of sim- 
ulated annealing [25[ and by pre-fixing to Landau gauge 
prior to MCG fixing [H. As observed in SU{2) [HUlf , 
there appears to be a significant anti-correlation between 
the value achieved in the gauge fixing functional and the 
percentage string tension reproduced by center vortices. 

As can be seen in Fig. [4] the ratio between the vortex- 
only string tension and the vortex density (simply the 
fraction of vortex plaquettes to total number of plaque- 
ttes) as a function of preconditioning is not independent 
of the preconditioning. Had it been independent one 
might conclude the reduction in the string tension is asso- 
ciated with simply not identifying all the vortices present. 
Either the mechanism with which vortices produce con- 
finement is not entirely intact or the physical relevance 
of the vortices is not uniformly distributed. 

In SU(2), it was seen that smearing an SU(2) configu- 
ration prior to MCG fixing reduced the center projected 
string tension considerably [1| . It was argued there that 
this is because smearing greatly expands the vortex cores 
making the MCG process of collapsing them to pierce a 
single plaquette more difficult. A similar point was used 
to address the issue raised by prefixing to Landau gauge 
(24j. In principle the same position could be taken here, 
the generated preconditioning transformation may allow 
the vortex cores to be distributed across a larger number 
of lattice sites and again make the MCG task of com- 
pressing them more difficult. However, the over-improved 
stout-link smearing parameters are deliberately chosen to 
maintain the size of instantons and there is a case to say 
that if the link between center vortices and topology seen 
in SU(2) persists in SU(3) then it should be possible to 
smear configurations without expanding the vortex cores. 
It is difficult to attribute the same vortex-expanding be- 
havior to the case of simulated annealing. However, the 
fact that known higher maxima exist (having these prop- 
erties) and that simulated annealing is designed to locate 
them could explain the similar behavior. 

It is significant that much of the discussion in Ref . [l5[ , 
where there is no Gribov ambiguity, can also be recon- 
ciled with the results found here. In this case a number 
of different Laplacian operators were constructed simply 
by using smeared links in the definition of the operator. 
There too it was seen that this caused an analogous ef- 
fect on the vortex-only string tension. It was argued that 
the use of smearing caused the Laplacian to be blind to 
the short-range physics making the decomposition of the 
gauge field into the confining and non-confining compo- 
nents less effective — disorder in the vortex-only com- 
ponent is absorbed into the vortex-removed component 
resulting in a loss of string tension. It was contended 
there that in the continuum limit the smearing radius 
shrinks to zero, restoring the string tension. 

In the same way, the smearing preconditioning may al- 
low this effect to occur for MCG. That the locations of 
vortices as determined by both methods coincide serves 
to strengthen this position. Indeed, as discussed in 
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Ref. [ll|, periodic boundary conditions cause gauge de- 
fects to have an opposite partner and, perhaps, the non- 
locality introduced by the preconditioning procedure may 
allow these opposites to annihilate producing no net de- 
fect after projection and a resultant drop in the vortex- 
only string tension. 

IV. CONCLUSIONS 

The use of smearing as a preconditioning technique 
leads to higher maxima in the MCG gauge-fixing condi- 
tion Vj/[f2]. These higher maxima in turn lead to lower 
numbers of P-vortices determined in the center projec- 
tion and, subsequently, lower values of the vortex-only 
string tension. 

Although the fundamental modular region of MCG 
would be an ideal candidate for a unique definition of 



vortex texture, it seems that the vortex matter arising 
from the first Gribov region as a whole has a greater 
phenomenological relevance. While all preconditionings 
lead to a loss of string tension, it is the center-projected 
physics that is not consistent. An improvement in Vj/[f2] 
causes one to miss vortices in the projection and spoil 
the phenomenology. 

While MCG has proved successful to a large extent 
in SU(2), what is different in 517(3) is that center- 
projection has never been shown to find enough vortices 
to reproduce the full string tension [38J. It just may 
be that MCG gauge fixing criterion is not sufficient to 
accurately locate the center vortices in SU(3) and that 
reproducing the success in SU(2) is a matter of getting 
closer to an "ideal" gauge. Nevertheless, it would be in- 
formative to look for correlations between the locations 
of the determined P-vortices with each preconditioning 
since their removal still leads to a loss of string tension. 
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